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Abstract 

In this article, we assume that the narrow structure X(4350) is a scalar 
cc — D*D* mixing state, and study its mass using the QCD sum rules. The 
i_j ' numerical result Mx = (4.37 ± 0.15) GeV is in good agreement with the 

j^. ■ experimental data, the X(4350) may be a scalar cc—D*D* mixing state. Other 

t-h , possibility, such as a scalar (tensor) cscs tetraquark state is not excluded. 

S PACS number: 12.39.Mk, 12.38.Lg 
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' — '! 1 Introduction 

In 2009, the CDF collaboration observed a narrow structure (which is denoted as the 
y(4140) now) near the J/ip(j) threshold with statistical significance in excess of 3.8 
standard deviations in exclusive decays B + — > J/ip(f)K + produced in pp collisions 
at y/s = 1.96 TeV pQ. The mass and width are (4143.0 ± 2.9 ± 1.2) Me V and 
(11.7jJ'o ± 3.7) MeV respectively. The narrow structure V(4140) is very similar to 
the charmonium-like state V(3930) near the J/ipoo threshold [2j [3]. The mass and 
width of the Y (3930) are (3914.6±|j ± 2.0) MeV and (34±* 2 ± 5) MeV respectively 



There have been several explanations for the nature of the narrow structure 
V(4140), such as a D*D* molecular state 01 El El El El ED] , an exotic (J PC = 1~+) 
hybrid charmonium [5], a cess tetraquark state |llj . the effect of the J/ip4> threshold 
[12] , or none a conventional charmonium state [32] nor a scalar D*D* molecular state 
[Til [15], etc. Assuming the V(4140) is a D*D* molecular state with J PC = ++ or 
2 ++ , Branz et al predict its two-photon decay width is of order 1 KeV [6j. 

Recently, the Belle collaboration measured the process 77 — > (pJ/ip f° r the 4>J/il> 
invariant mass distributions between the threshold and 5 GeV based on a data sam- 
ple of 825fb~ 1 , and observed a narrow peak of 8.8I32 events with a significance 
of 3.2 standard deviations [TB]. The mass and width of the structure (denoted as 
X(4350)) are (4350.6±|;f ± 0.7) MeV and (13.3^ 7 £ 9 ± 4.1) MeV respectively. No sig- 
nal for the V(4140) — > <pJ/ip structure was observed, this disfavors the scenario of 
the V(4140) as a D*D* molecular state. 

The possible quantum numbers for a state X decaying into J/ip<j) are 

jPC = 

0~ + , ++ , 1 _+ , 2 ++ , the corresponding strong interactions can be described by the 
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following phenomenological Lagrangian, 




(1) 



the strong coupling constants g can be fitted phenomenologically or calculated by 
some theoretical approaches, for example, the QCD sum rules. 

In Ref. [5], Zhang and Huang study the QsQs and QsQ's molecular states in a 
systematic way using the QCD sum rules before the Belle experiment, the mass of the 
D* + Dj molecular state is (4.36 ± 0.08) GeV, which is consistent with experimental 
data [IB] . Such a state has J p = 1~ and no definite charge conjugation. In Ref . [IT] . 
Albuquerque et al re-study the exotic D* + D~ — D*~Df molecular state with J PC = 
1~ + by taking into account the contributions from the vacuum condensates up to 
dimension-8, the prediction Md*d s0 = (5.05 ± 0.19) GeV disfavors the scenario of 
the X(4350) as a D* s D s q molecular state. 

In Ref. [18] , Drenska et al study the exotic tetraquark states of the kind cscs by 
computing their spectrum and decay modes within a constituent diquark-antidiquark 
model, the predictions M -+ = 4277, 4312 MeV and M x -+ = 4321, 4356 MeV are 
consistent with the experimental data [TE] . On the other hand, the flux-tube model 
[T9l [20] and the Lattice QCD j2H I23 123] predict that the masses of the low lying 
hybrid charmonium states are about (4.0—4.2) GeV and (4.0—4.4) GeV respectively, 
which are also consistent with the experimental data [16] . However, the decay of a 
hybrid to two photons is generically forbidden [24] . 

In Ref. [TTj . Stancu study the mass spectrum of the cess tetraquarks using a 
simple quark model with chromo-magnetic interaction and observe that the V (4140) 
may be the strange partner of the V(3872), the prediction for the mass of the 
2 ++ tetraquark state is consistent with the V(4350). As noticed by the author, 
the amplitude of the singlet-singlet component seems too large comparing with the 
octet-octet component. 

In Ref. [25] . Liu et al discuss the possibility that the X(4350) is an excited P- 
wave charmonium state Xc2 by studying the strong decays of the P-wave charmonium 
states with the 3 Po model. 

The CDF and Belle collaborations analyze the experimental data by assuming 
the vector mesons J/ip and have a relative S"-wave [TJ [T6J, so we will not focus 
on the scenarios of the V(4350) as the 0~ + and 1 h tetraquark state or hybrid 
charmonium. 

In Refs. [261 27], we study the mass spectrum of the scalar hidden charm and 
hidden bottom tetraquark states which consist of the scalar-scalar type, axial- vector- 
axial- vector type and vector- vector type diquark pairs in a systematic way using the 
QCD sum rules, the scalar-scalar type and axial-vector-axial-vector type scalar cess 
tetraquark states have masses about (4.45 ± 0.16) GeV, the lower bound of the 
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Table 1: The possible explanations for the nature of the X(4350), where the RW 
stands for the relative wave of the final state mesons. 



masses is consistent with the X(4350), we cannot exclude that the X(4350) is a 
scalar cess tetraquark state. In Refs.[T41 [T5]. we study the D*D*, D*D*, B*B* and 
B*B* molecular states in a systematic way using the QCD sum rules. The numerical 
result is inconsistent with the experimental data, the D*D* is probably a virtual 
state and not related to the meson F(4140). In this article, we study the X(4350) 
as a linear superposition of a scalar charmonium state cc and a virtual state D*D* 
using the QCD sum rules [2HI ES] • In Table 1, we present the possible explanations 
for the nature of the X(4350). 

The article is arranged as follows: we derive the QCD sum rules for the nar- 
row structure X(4350) in Sect. 2; in Sect. 3, we present the numerical results and 
discussions; and Sect. 4 is reserved for our conclusions. 



2 QCD sum rules for the X(4350) as a mixing state 

In the following, we write down the two-point correlation function H(p) in the QCD 
sum rules, 

n(p) = i J d 4 xe ip x (0\T{j(x)j\0)} |0), (2) 
_ Ji(x) + J 2 (x) 

Ji(x) = c(x)7 A1 s(x)s(x)7 A 'c(x) , 

J2(x) = ———c(x)c(x) , (3) 



where the J2(x) is the normalized two-quark current [30J. 

We can insert a complete set of intermediate hadronic states with the same 
quantum numbers as the current operator J(x) into the correlation function Tl(p) 
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to obtain the hadronic representation [281 EH]. After isolating the ground state 
contribution from the pole term of the lowest state X, we get the following result, 



\ 2 



M\-p 2 



+ ■ 



where the pole residue (or coupling) is defined by 

X x = (0|J(0)|X(p)>. 



(4) 



(5) 



The two-quark current c(x)c(x) has non- vanishing coupling with the charmonia 
X c o(lP), Xco(2-P), Xco(3-P), • • • ; while the molecule type current c(x)j IJi s(x)s(x) r y ,J 'c(x) 
has non- vanishing coupling with the molecular states D*D*, D*D'* S , D'* S D' S , ■ ■ ■ and 
the scattering states D* — D*, D* — D'* s , D'* s — D'* s , ■ ■ ■ [31] • We cannot distinguish 
those contributions and study them exclusively. In this article, we take the assump- 
tion that the interpolating current J(x) couples to a particular resonance, which is a 
special superposition of the scalar charmonia Xco(l-P); • • • and the virtual molecular 
states D*D*, ■ ■ ■ . In other words, we take a single pole approximation, the pole 
embodies the net effects. 

We carry out the operator product expansion for the correlation function Tl(p) 
at the large space-like momentum region p 2 CO, 



1 (ss) 2 
n(p) = -U n (p) + { -^U 22 (p) 



(6) 



where 
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n 22 (p) 
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7T 



dx 
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x' 



— x + l)ml 
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(7) 



the explicit expression of the spectral density pn(s) can be found in Refs.fTU [15], 

A 2 = 4(m c +m s ) 2 , fh 2 



I + a/1- 



4m 2 



/2, Xi 



1 - a/1 - 



4m? 



/2- 



x{X-x)i X f 

In calculation, we use the Fierz re-ordering in the color space and Dirac spin space 
to express the correlation functions Ilulp) and n. 2 i(p) in terms of the H22G9). In 
this article, we carry out the operator product expansion to the vacuum condensates 
adding up to dimension-10 and take the assumption of vacuum saturation for the 
high dimensional vacuum condensates, they are always factorized to lower conden- 
sates with vacuum saturation in the QCD sum rules, and factorization works well 
in large iV c limit. 
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Once analytical result is obtained, then we can take the quark-hadron duality 
and perform the Borel transform with respect to the variable P 2 = —p 2 , finally we 
obtain the following sum rule: 



ds 



A2 



2P110) + 



>ssy 

6 



-P22{S) 



e 



(8) 



Differentiating Eq.(8) with respect to j^, then eliminate the pole residue \x, 
we can obtain the sum rule for the mass, 



Ja° ds d(-i/Mi) [ 3 Pn( g ) + ( ss ) 2 P22(s)} e ^ 
J^ 2 ds [3pu(s) + (ss) 2 p 2 2(s)} e"i^ 



(9) 



3 Numerical results and discussions 

The input parameters are taken to be the standard values (qq) = — (0.24±0.01 GeV) 3 , 
(ss) = (0.8 ± 0.2) (qq), (sg s aGs) = m 2 (ss) } m 2 = (0.8 ± 0.2) GeV 2 , = 
(0.33 GeV) 4 , m s = (0.14 ± 0.01) GeV and m c = (1.35 ± 0.10) GeV at the energy 
scale /1 = 1 GeV [281 [291 [32] . 

In the conventional QCD sum rules [28j [29], there are two criteria (pole dom- 
inance and convergence of the operator product expansion) for choosing the Borel 
parameter M 2 and threshold parameter so- In Refs. |141 IT5] . we study the D*D*, 
D*D*, B*B* and B*B* molecular states in a systematic way, the threshold parame- 
ters are s = (24 ±1) GeV 2 , (25 ±1) GeV 2 , (138 ±2) GeV 2 and (140 ±2) GeV 2 in the 
c r y fl ud , j fl c, c r y^ss'y fl c, b^^ud^b and b^^ss^b channels respectively; the Borel param- 
eters are M 2 = (2.6-3.0) GeV 2 and (7.0-8.0) GeV 2 in the hidden charm and hidden 
bottom channels respectively. In those regions, the two criteria of the QCD sum 
rules are satisfied. In this article, we choose the interpolating current J(x), which is 
a special superposition of the scalar currents c(x)c(x) and c(x)7 /i s(x)s(x)7 M c(x). So 
we can take the same threshold parameter and Borel parameter as in the channel 
c^ss^c, i.e. s = (25 ± 1) GeV 2 and M 2 = (2.6 - 3.0) GeV 2 . 

The contributions from the different terms in the operator product expansion are 
shown in Fig.l, from the figure, we can see that the dominant contribution comes 
from the perturbative term and the operator product expansion is well convergent. 
In Fig.2, we show the contribution from the pole term with variation of the Borel 
parameter and the threshold parameter. The pole contribution is larger than 50%, 
the pole dominant condition is also satisfied. 

Taking into account all uncertainties of the relevant parameters, finally we obtain 
the values of the mass and pole residue of the narrow structure V(4350), which are 
shown in Fig. 3, 

M x = (4.37 ±0.15) GeV, 

X x = (4.1 ±0.8) x 10~ 2 GeV 5 . (10) 
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The prediction is in good agreement with the experimental data Mx = (4350.6l 51 ± 
0.7) MeV [To] , the X(4350) may be a scalar cc — D*D* mixing state. Other possi- 
bility, such as a scalar (tensor) cess tetraquark state is not excluded. 

The nominal thresholds of the D s — D s and D* — D* are M Da [) s = 3.937 GeV and 
M D , D , = 4.225 GeV respectively [3T], the strong decays X(4350) ->• D S D S ,D*D* 
can take place, so we can search for the X(4350) in the D S D S and D*D* invari- 
ant mass distributions. The decay channel X(4350) — > D*D* has much smaller 
phase space comparing with the decay channel X(4350) — > D S D S , the strong decay 
X(4350) — > D S D S is of great importance. By measuring the relative angular distri- 
butions of the pseudoscalar mesons D s and D s , we can determine the spin of the 
X(4350). 

In Ref . [33] , the light nonet scalar mesons are taken as tetraquark states, and the 
strong coupling constants among the light scalar mesons and pseudoscalar mesons 
are calculated with the QCD sum rules. The numerical results indicate that the 
values of the strong coupling constants for the tetraquark states are always smaller 
than the corresponding ones for the qq states [MJ [35]. In Ref. [36], Maiani et al take 
the diquarks as the basic constituents, examine the rich spectrum of the diquark- 
antidiquark states with the constituent diquark masses and the spin-spin interac- 
tions, and try to accommodate some of the newly observed charmonium-like reso- 
nances not fitting a pure cc assignment. The predictions (also the Ref. [18]) depend 
on the assumption that the light scalar mesons a (980) and /o(980) are tetraquark 
states, the basic parameters (constituent diquark masses) are estimated thereafter. 
If the scenarios of the light nonet scalar mesons as the tetraquark states are robust, 
the scalar (tensor) cess tetraquark state will have smaller D S D S decay width than 
the corresponding ones of the cc — D*D* mixing state. 

4 Conclusion 

In this article, we assume that the X(4350) is a scalar cc — D*D* mixing state, and 
study its mass using the QCD sum rules. Our prediction depends heavily on the two 
criteria (pole dominance and convergence of the operator product expansion) of the 
QCD sum rules. The numerical result is in good agreement with the experimental 
data, the X(4350) may be a scalar cc — D*D* mixing state. Other possibility, 
such as a scalar (tensor) cess tetraquark state is not excluded. We can search for 
X(4350) in the D S D S and D*D* invariant mass distributions, especially the D S D S . 
By measuring the relative angular distributions of the pseudoscalar mesons D s and 
D s , we can determine the spin of the X(4350). 
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Figure 1: The contributions from the different terms with variation of the Borel 
parameter M 2 for so = 25 GeV 2 in the operator product expansion. The A, 
B, C, D, E and F correspond to the contributions from the perturbative term, 
(ss) + (sg s aGs) term, ( s ^) term, ( Bm ^)+( 9 ^) [(as) + (sg s aGs) + (ss) 2 } term, 
(ss) 2 +(ss)(sg s aGs) term and (sg s aGs) 2 term, respectively. 
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